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ABSTRACT 

We correlate the positions of radio galaxies in the FIRST survey with the CMB lens- 
ing convergence estimated from the Atacama Cosmology Telescope over 470 deg^ to 
determine the bias of these galaxies. We remove optically cross-matched sources below 
redshift 2 = 0.2 to preferentially select Active Galactic Nuclei (AGN). We measure 
the angular cross-power spectrum at 4.4(j significance in the multipole range 
100 < I < 3000, corresponding to physical scales between 2-60 Mpc at an effective 
redshift Zeff = 1 - 5 . Modelling the AGN population with a redshift-dependent bias, the 
cross-spectrum is well fit by the Planck best-fit AGDM cosmological model. Fixing the 
cosmology we fit for the overall bias model normalization, finding 6(2eff) = 3.5zb0.8 for 
the full galaxy sample, and b{zeif) = 4.0 zb 1.1 (3.0 zb 1.1) for sources brighter (fainter) 
than 2.5 mJy. This measurement characterizes the typical halo mass of radio-loud 
AGN: we find log(Mhaio/Af0) = 13 . 6 ^ 0 ^. 

Key words: large-scale structure of Universe, cosmic microwave background, radio 
continuum: galaxies. 


1 INTRODUCTION 

Radio galaxies trace the large-scale structure in the Universe 
which has been measured with large-area surveys includ- 
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ing FIRST, WENSS, NVSS, and SUMSS (Becker, White & 
Helfand 1995; Rengelink et al. 1997; Condon et al. 1998; 
Bock, Large & Sadler 1999); for an overview see de Zotti 
et al. (2010). The angular clustering of these galaxies has 
been measured by Cress et al. (1996); Magliocchetti et al. 
(1998); Blake & Wall (2002); Overzier et al. (2003); Blake, 
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Mauch & Sadler (2004); Lindsay et al. (2014). The cluster¬ 
ing of radio galaxies will soon be measured over much larger 
volumes of the universe with the Square Kilometer Array 
(SKA) and its precursors, allowing cosmological effects such 
as dark energy, modified gravity and non-Gaussianity to be 
probed (e.g., Carilli & Rawlings 2004; Blake et al. 2004; Rac- 
canelli et al. 2012; Camera et al. 2012; Maartens et al. 2013; 
Norris et al. 2013; Jarvis et al. 2015). 

The bias 6 of a large-scale structure tracer relates over¬ 
densities of that tracer S to overdensities of the underlying 
dark matter field Jdm: 

5 — bSoM- (1) 

Radio-selected galaxies broadly contain two populations: 
high-redshift Active Galactic Nuclei (AGN) and low-redshift 
star-forming galaxies (Gondon, Gotton & Broderick 2002). 
AGN dominate the radio emission at high flux (> 1 mjy) 
and are highly biased, their hosts being among the most 
massive galaxies in the early universe (e.g., Jarvis et al. 
2001a; Rocca-Volmerange et al. 2004; Seymour et al. 2007; 
de Zotti et al. 2010; Fernandes et al. 2015). Their bias de¬ 
pends strongly on galaxy mass and redshift (e.g., Seljak & 
Warren 2004), and is poorly constrained particularly at high 
redshift where few optical counterparts are observed. Some 
progress has been made by identifying redshifts spectroscop¬ 
ically: using Galaxy And Mass Assembly (GAMA) data the 
bias of FIRST radio galaxies was measured at « « 0.34 over 
200 deg^ to the 10% level (Lindsay et al. 2014). On a smaller 
square degree region, clustering measurements using data 
from the Very Large Array (VLA) and VISTA Deep Ex- 
tragalactic Observations (VIDEO; Jarvis et al. 2013) were 
used to show evidence for a strongly increasing bias a.t z > 2 
(Lindsay, Jarvis & McAlpine 2014). 

An alternative way to constrain bias is through cross¬ 
correlation of the tracer fluctuations with gravitational leas¬ 
ing due to large-scale structure. In particular, the leasing of 
the Cosmic Microwave Background (CMB) measures the in¬ 
tegrated matter fluctuations to 2 ; « 1100. As we will show 
the high-redshift radio source distribution overlaps strongly 
with the broad CMB leasing kernel. Cross-correlations be¬ 
tween the CMB and other tracers of large-scale structure 
have been reported by e.g.. Smith, Zahn & Dore (2007); Hi- 
rata et al. (2008); Feng et al. (2012); Bleem et al. (2012); 
Planck Collaboration (2014c); van Engelen et al. (2014); 
Fornengo et al. (2014). The Planck Collaboration (2014b) 
detect the correlation of leasing with radio galaxies from 
NVSS at 20(7. Sherwin et al. (2012) correlate leasing mea¬ 
surements from the Atacama Cosmology Telescope (ACT) 
with optically-selected quasars from the Sloan Digital Sky 
Survey (SDSS), measuring a bias of b = 2.5 ± 0.6 at an ef¬ 
fective redshift « « 1.4. Geach et al. (2013) correlate lensing 
from the South Pole Telescope (SPT) with quasars selected 
from the Wide-held Infrared Survey Explorer (WISE), mea¬ 
suring a bias b = 1.61 ± 0.22 at 2 « 1.0. An advantage 
of cross-correlations is that they are robust to systematic 
biases which may be particular to each dataset. 

In this paper we measure the angular cross-power spec¬ 
trum between the lensing convergence estimated from 
AGT with the FIRST radio source overdensity. We use lens¬ 
ing maps from the three-year ACT Equatorial survey (Das 
et al. 2013) together with the hrst-season ACTPol survey 
(Naess et al. 2014; Madhavacheril et al. 2014; van Enge¬ 


len et al. 2014). We consider 36,000 radio sources with flux 
brighter than 1 mJy, and remove optically cross-matched 
sources from the Sloan Digital Sky Survey (York et al. 2000) 
at 2 < 0.2 to preferentially select AGN, discarding the ma¬ 
jority of low-redshift star-forming galaxies. We use this to es¬ 
timate the bias normalization, assuming a hxed cosmological 
model, and using a redshift distribution and bias-evolution 
model from the simulated radio catalogue of the SKA De¬ 
sign Study (SKADS, Wilman et al. 2008). We measure Gf® 
across a wide range of scales (100 < I < 3000) and consider 
various splits of the radio sources to investigate redshift and 
flux dependence of the bias. 

We describe the lensing and radio data and the cross¬ 
correlation analysis methods in Section 2. The results and 
discussion are presented in Section 3, with further inter¬ 
pretation of the AGN bias in Section 3.1. We conclude in 
Section 4. 

2 DATA AND ANALYSIS 
2.1 ACT and ACTPol 

The Atacama Cosmology Telescope (ACT) is located at 
an altitude of 5190m in Parque Astronomico Atacama in 
Northern Chile. The telescope and its current polarization- 
sensitive receiver, ACTPol, are described in Niemack et al. 
(2010). The two seasons of ACT temperature data and the 
ACTPol first-season temperature and polarization data used 
in this analysis are presented in Das et al. (2013) and Naess 
et al. (2014). Lensing by large-scale structure induces cou¬ 
pling of otherwise independent temperature and polarization 
modes. We construct estimators of the lensing convergence 
from quadratic combinations of temperature and polariza¬ 
tion maps in Fourier space, following the methodology of 
Hu & Okamoto (2002). We use the same lensing conver¬ 
gence maps and Monte-Carlo simulations described by Das 
et al. (2011) and van Engelen et al. (2014). 

In this analysis we use two ACT datasets. The first is 
the ACT Equatorial data which spans a thin strip along 
the celestial equator with an area of 300 deg^. This strip is 
partitioned into six approximately equal area patches over 
which we compute the cross-spectrum separately then aver¬ 
age (weighting by patch area) for the final result. We lose 
negligible information at the scales of interest and it allows 
for patch to patch consistency checks. The effective white- 
noise component of the two-season co-added data is 18 jzK- 
arcmin. 

We also fold in the three ACTPol ‘deep’ fields from 
the first-season dataset, labelled Dl, D5 and D6, with a 
temperature white-noise component of 16.2,13.2,11.2 /rK- 
arcmin, respectively, over a total area of 206 deg^ (37 deg^ 
of which overlaps with the ACT Equatorial strip). All maps 
in this analysis use 0.5' x 0.5' pixels, and their positions are 
shown in Fig. 1. 

For each ACTPol patch we use the minimum-variance 
(MV) linear combination of the reconstructed convergence 
maps estimated from each quadratic pair (TT, TE, EE, EB). 
Following the POLARBEAR Collaboration (2013), van En¬ 
gelen et al. (2014) and Story et al. (2014) we perform the 
combination in Fourier space, weighting each convergence 
map by the mode-dependent inverse-variance noise to ob¬ 
tain the MV combination. All lensing convergence maps are 
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mean-field subtracted to remove the lensing-like effect at 
large scales of mode-coupling from the windowing of the 
temperature and polarization fields. 

As described in van Engelen et al. (2014) an apodization 
window is applied to the ACT and ACTPol temperature and 
polarization map prior to lensing reconstruction. This win¬ 
dowing operation includes a cosine taper at the map edges 
to remove discontinuous edges and weighting by the pixel 
bitmap to optimize the signal-to-noise of the reconstruction. 
The resulting quadratic estimator reconstruction is therefore 
also windowed, resulting in a scale-dependent suppression of 
power. Following Bleem et al. (2012), Sherwin et al. (2012), 
Hand et al. (2013) and van Engelen et al. (2014), we use 
realistic Monte-Carlo simulations to calculate the transfer 
function correction by computing the mean cross-spectrum 
between noiseless lensing realizations and their correspond¬ 
ing reconstructions within the lensing pipeline. This correc¬ 
tion (< 5% for ACTPol, « 10% for ACT) is then applied 
to the maps when computing the data cross-spectrum to 
account for the suppression of power due to windowing. 

2.2 FIRST 

The FIRST survey (Becker, White & Helfand 1995) was car¬ 
ried out between 1993 and 2011 at 1.4 GHz with the VLA 
in B configuration. The final catalogue (Helfand, White & 
Becker 2015) contains 946,432 sources covering 10,575 deg^, 
with an angular resolution of 5.4” (FWHM) and to a com¬ 
pleteness of 95% at flux Si, 4 GHz > 2 mjy. The oblique 
decision-tree program developed by the FIRST survey team 
(White et al. 1997) determines the probability that each cat¬ 
alogue entry is the result of a spurious sidelobe response to 
a nearby bright source. We exclude entries with a sidelobe 
probability of > 0.1, leaving 720,219 sources above 1 mJy. 

To address the issue of extended radio sources resulting 
in multiple detections for one host galaxy, perhaps none of 
which corresponds to the core itself (and therefore any asso¬ 
ciated optical source), we have followed Cress et al. (1996) in 
applying a collapsing radius of 72” (0.02 deg) to the FIRST 
catalogue. Any FIRST sources within this radius of one an¬ 
other are grouped and combined to form a single entry, posi¬ 
tioned at the flux-weighted average coordinates of the group, 
and attributed with their total flux density. Around 32% of 
all FIRST sources are collapsed (in groups of average size 

2.3 sources per group), forming 17% of the resulting cat¬ 
alogue. These multiple-component sources will come from 
AGN, which dominate the source population at high flux 
density and high redshift; at lower flux density and redshift, 
starbursts and normal star-forming galaxies (SFGs) are in¬ 
creasingly dominant (Condon, Cotton & Broderick 2002). 

The redshifts of individual sources are not determined 
by FIRST, but can be found by identifying counterparts in 
SDSS, which gives redshifts for the brighter, nearby sources 
in the FIRST sample. The closest sources are most likely 
star-forming galaxies (e.g., Condon, Cotton & Broderick 
2002; Wilman et al. 2008); by removing them we simplify 
the measurement as a constraint on the bias of the domi¬ 
nant astrophysical population (i.e. AGN). 

We do this by initially taking all sources in the cata¬ 
logue which lie in the ACT and ACTPol patches described 
above (« 38,000 sources). AGN dominate the radio luminos¬ 
ity function for Li, 4 ghz > 10^^ W Hz~^ (Condon, Cotton 


& Broderick 2002; Jarvis & Rawlings 2004; Mauch & Sadler 
2007). Given the flux limit of FIRST (ImJy), and assuming 
a spectral index of a = 0.8 {Si, oc for the AGN, this 
luminosity threshold corresponds to sources above redshift 
« = 0.2. We identify optical matches to the radio sources 
within SDSS, treating as reliable all matches within 2” of the 
radio source following Lindsay et al. (2014). Given the den¬ 
sity of SDSS sources, the level of spurious optical matches 
identified with this technique is below 2%. A fraction of 0.27 
of the FIRST sources in the lensing fields have an optical 
match obtained in this manner. We remove all sources with 
a known redshift below « = 0.2, constituting 18% of the 
sources with a reliable redshift, or 5% of the total number 
of sources. Given the small fraction of sources removed, this 
procedure has only a small effect on the results (Section 3). 
The final sample comprises « 36,000 sources with a mean 
angular density of 71 sources deg“^. 

Within each ACT and ACTPol patch a corresponding 
map of the overdensity of sources g is produced in a similar 
way to Sherwin et al. (2012) and Geach et al. (2013). We 
create a map at the same pixelation as the lensing map and 
define the radio-galaxy overdensity map g by 


where rii is the number of sources in each pixel and h is the 
mean number of galaxies per pixel. This overdensity map is 
then smoothed with a Gaussian with FWHM of 2' to obtain 
a well-defined pixel window function. 

2.3 Analysis methods 

We compute the cross-spectrum between the lensing conver¬ 
gence from AGT and ACTPol with the FIRST radio galaxy 
overdensity. 

Following the procedures outlined in Das et al. (2011) 
and Hand et al. (2013), we correct for mode-coupling in¬ 
duced by windowing in real space and from applying annular 
binning in Fourier space, computing an unbiased estimator 
of the binned cross-spectrum The binning we adopt is 
given in Table 1. 

To determine the full band-power covariance matrix 
we cross-correlate realistic simulations of the reconstructed 
lensing fields with the radio source maps (which are in prin¬ 
ciple uncorrelated). Production of these realistic simulations 
is described in Das et al. (2011) and van Engelen et al. 
(2014). 

This procedure ignores the cosmic variance contribution 
to the uncertainties in the data coming from the correlated 
part of the two maps, C'C ■ We neglect this as both maps are 
noise-dominated at the relevant scales for this analysis. Bin- 
to-bin correlations are < 10% for all off-diagonal elements of 
the covariance matrix. We also check that the mean cross¬ 
spectrum is consistent with null (Figure 2), confirming that 
our pipeline does not induce spurious cross-power in the 
absence of correlation. 

Approximately 50% of the ACTPol D5 patch and 15% 
of D6 overlap with the ACT Equatorial strip (Fig. 1). There 
is therefore a correlation between the ACT and ACTPol 
cross-spectra, as common CMB modes in the primary tem¬ 
perature map have been used to reconstruct the lensing 
convergence over these regions. Noiseless temperature maps 
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Figure 1. Top panel: Footprint of the patches used in this anal¬ 
ysis: ACT (blue) and ACTPol (green, left to right: Dl, D6, D5). 
Middle panel: ACTPol Dl lensing convergence map k smoothed 
to suppress power below 20' scales. The spatial modulation, pri¬ 
marily due to the windowing of the temperature and polarization 
maps by the pixel weight map, is evident. The ACTPol lens¬ 
ing convergence is noise dominated for scales < 1 degree. Bot¬ 
tom panel: The FIRST overdensity field g over the same patch, 
smoothed to the same scale, is noise dominated at all scales. 


from ACT and ACTPol, and negligible polarization infor¬ 
mation from ACTPol, would result in a perfect correlation 
between the reconstructed convergence maps. However, this 
overlapping area represents 37 deg^ of the total 470 deg^ of 
this analysis, and hence at most a 4% overestimate of the 
detection significance, which we neglect given the statistical 
errors. We thus average the ACT and ACTPol data cross¬ 
spectra with inverse-variance weighting. 

In order to check for bias in the cross-spectrum estima¬ 
tor, we ran 500 pairs of simple simulated convergence and 
radio density maps through the cross-correlation pipeline, 
generating new correlated simulations. To obtain these pairs 
we draw as signal maps aperiodic correlated Gaussian real¬ 
izations from power spectra obtained assuming Planck best- 
ht cosmological parameters and a fiducial bias model and 
source distribution for the radio galaxies (Kamionkowski, 
Kosowsky & Stebbins 1997; Planck Collaboration 2014a). 
We add Gaussian noise realizations to the convergence maps, 
appropriate for the temperature sensitivity of ACT (Section 
2.1), using the formalism of Hu & Okamoto (2002) to calcu¬ 
late the reconstruction noise. ACTPol maps are less noisy, 
but the precise noise level is unimportant for this test. For 
each pixel i in the radio signal map g we draw a Poisson ran- 



Figure 2. Pipeline validation: Mean cross-spectrum IC^^ of the 
FIRST radio source map with ACTPol lensing simulations (green 
dashes, Aigims = 2048) and ACT lensing simulations (blue circles, 
A'sims = 480) as described in 2.3. We displace the ACT points 
by Al = 30 to the right for visual clarity. The measurements are 
consistent with null, demonstrating that our pipeline does not in¬ 
duce spurious cross-power in the absence of correlation. Error bars 
shown are the diagonal components of the empirical covariance 
matrix derived from the same Monte Carlo simulations, scaled 
appropriately by We also show the recovered mean cross¬ 

spectrum from realistic correlated simulations (red triangles. Sec¬ 
tion 2.3). We cross-correlate input convergence maps, which have 
added scale-dependent Gaussian noise, with correlated realiza¬ 
tions of a galaxy field. This demonstrates that our pipeline is able 
to recover in an unbiased fashion a known input cross-spectrum 
(although we note this does not test the lensing reconstruction 
pipeline, for which we refer to the systematic tests in van Enge- 
len et al. (2014)). The generative model for the cross-spectrum is 
not the fiducial cross-spectrum, but this is unimportant for the 
purposes of this test. 


dom variable Xi with mean n(l-\-gi), where n is the average 
number of sources per pixel. We set h = 71 sources deg“^ 
to reflect the source density in the data. We then redefine 
gii - 1 — Xi/h — 1 and Hnally smooth the resulting map with a 
Gaussian beam of FWHM 2'. 

These simulated maps, by construction, have signal, 
noise and correlation properties which mimic the data, al¬ 
though they do not have the full spatially anisotropic noise 
properties. These lensing simulations have not been pro¬ 
cessed through the lensing reconstruction pipeline, but here 
we use them simply for checking bias in the cross-correlation 
pipeline. We refer to the systematic tests in van Engelen 
et al. (2014) for checks of the lensing pipeline. We find that 
we do not require apodization of the maps to produce the 
observed unbiased results; the mean auto- and cross-spectra 
of these simulations are consistent with the assumed input 
spectra (Fig. 2). 


2.4 Modelling 

The theoretical cross-spectrum can be written under the 
Limber approximation as 

where H(z) is the Hubble parameter, x(A i® ^he comov¬ 
ing distance to redshift z, P{k, z) is the non-linear matter 
power spectrum (wavenumber k — l/x) and {Wi} are the 
appropriate kernels for the two dark matter probes k, g. The 
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Figure 3. Green dashed: Fiducial bias model b{z), constructed as 
a weighted average of the constituent source population bias mod¬ 
els (Section 2.2). Blue solid: Source redshift distribution model as 
derived from SKADS and including the cut of a fraction of the 
2 < 0.2 sources. Vertical black line: z = 0.2. Grey solid: As previ¬ 
ous but without the redshift cut. Red dot-dashed: Unnormalised 
CMB lensing kernel Wk{z). 


dominant term in Wg is directly proportional to the tracer 
bias b{z): 

WAz) = b{z)^+M{z), (4) 

where dn/dz is the normalised source redshift distribution 
and M{z) is a sub-dominant contribution from the mag- 
nihcation bias (see e.g., Sherwin et al. (2012) for the full 
expression). The magnification bias term is independent of 
the tracer bias, and for the FIRST sources in our sam¬ 
ple is small (< 6 % of the total). A rescaling of the bias 
amplitude therefore corresponds linearly to a rescaling of 
the cross-spectrum . We compute the theory P{k, z) us¬ 
ing best-fit Planck cosmological parameters, including non¬ 
linear corrections using CAMB with Halofit (Lewis, Challi- 
nor & Lasenby 2000; Smith et al. 2003; Takahashi et al. 
2012 ). 

We use the SKA Design Study (SKADS) simulated ra¬ 
dio continuum catalogue to construct a fiducial bias model 
b{z) and redshift distribution dn/dz for the 5 i, 4 ghz > 1 mjy 
radio sources (see Wilman et al. 2008, 2010, for details). 
The simulation lacks the mass resolution to directly resolve 
galaxy- and group-sized haloes for a robust implementation 
of the galaxy clustering, but the source counts, redshift dis¬ 
tribution, and variations in space density are defined by ex¬ 
trapolating observed luminosity functions, and implement¬ 
ing a bias model, for each of hve individual radio popu¬ 
lations: AGN (FRI and FRII types, radio-quiet quasars), 
normal star-forming galaxies and starburst galaxies. These 
populations are assigned a single halo mass each, used to de¬ 
fine b{z) as described by Mo & White (1996), with the bias 
held fixed above a particular redshift to prevent unphysical 
clustering where the bias is poorly constrained observation- 
ally (see Figure 3 of Raccanelli et al. 2012). The simulated 
catalogue informs us how the relative numbers of these pop¬ 
ulations evolve with redshift, and how the observed bias will 
evolve accordingly for a mixed sample of sources. 

By comparing the distribution of known source red- 
shifts (Section 2.2) with the SKADS simulation, we hnd an 
estimated 66 % of low-redshift sources are removed by the 


Bin b 

[^min ! 

^max] 

^K.g 

^b,ACT 

(xlO'^) 

^Kg 

ACTPol 

(xlO'^) 

CC K 

o,comb 

(xlO'^) 

200 

[100,300] 

1.76 ±0.74 

2.04 ±0.80 

1.89 ±0.54 

450 

[301,600] 

0.59 ±0.32 

-0.20 ±0.42 

0.30 ± 0.25 

750 

[601,900] 

0.57 ± 0.25 

0.32 ±0.31 

0.48 ±0.19 

1050 

[901,1200] 

0.28 ±0.20 

0.38 ±0.26 

0.32 ± 0.16 

1350 

[1201,1500] 

0.17 ± 0.16 

-0.08 ±0.22 

0.09 ±0.13 

1650 

[1501,1800] 

0.04 ± 0.15 

0.01 ± 0.19 

0.03 ±0.12 

1950 

[1801,2100] 

0.11 ± 0.14 

-0.00 ±0.18 

0.07 ±0.11 

2250 

[2101,2400] 

0.07 ± 0.14 

-0.01 ± 0.17 

0.04 ±0.11 

2550 

[2401,2700] 

0.06 ±0.14 

0.19 ±0.18 

0.11 ± 0.11 

2850 

[2701,2999] 

-0.01 ± 0.13 

-0.23 ±0.19 

-0.08 ±0.11 


Table 1. The measured cross-spectrum for FIRST radio 

sources with ACT and ACTPol lensing. The bins are chosen to 
be wide enough that correlations are small (< 10%, Section 2.3), 
but narrow enough to resolve structure in the cross-spectrum. 



Figure 4. Data cross-spectrum 1C”® for (ACT-1-ACTPol) X 
FIRST. Dark grey, solid: the best-fit cross-spectrum. Light grey, 
solid: the contribution from z > 1.5 sources. We restrict the plot 
to I < 2000 where the signal-to-noise dominates. We show as error 
bars the diagonal components of the empirical covariance matrix 
derived from Monte Carlo simulations (Section 2.3). Scaling the 
amplitude of the fiducial bias model, the combined significance of 
the bias detection is 4.5 (t (Section 3). 

2 < 0.2 cut. To construct dn/dz we therefore weight 2 < 0.2 
sources by 0.34 relative to higher-redshift sources, similar to 
the approach of Lindsay, Jarvis & McAlpine (2014). After 
the redshift cut « 96% of the sources in our sample are ex¬ 
pected to be AGN, with a « 4% contamination fraction of 
star-forming and starburst galaxies, and we estimate the h- 
nal sample to have a median redshift 2 = 1.3. These fiducial 
models are shown in Figure 3. We discuss the limitations 
of this model, including the effect of not removing the low- 
redshift sources, in Section 3.3. 

Finally, we bin the theoretical cross-spectrum as for the 
data, accounting for the mode-coupling matrix of each patch 
(Das et al. 2011). We then compare the model to the data 
using a simple Gaussian likelihood, and primarily fit for an 
overall scaling A to the fiducial bias model, such that 6 ( 2 ) —>■ 
Ab{z). 

3 RESULTS 

The cross-spectra for ACT x FIRST, ACTPol x FIRST, 
and their combination are shown in Figs. 4 and 5 and 
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Figure 5. Data cross-spectrum for ACT X FIRST (blue) 

and ACTPol X FIRST {green). Dark grey, solid: the best-fit cross¬ 
spectrum for the combined data. Light grey, solid: the contribu¬ 
tion from z > 1.5 sources. The ACT and ACTPol points have 
been displaced to the right and left by Al = 15, respectively, for 
visual clarity. We show as error bars the diagonal components of 
the empirical covariance matrix derived from Monte Carlo simu¬ 
lations (Section 2.3). 



A 

S/N 

(A 

PTE 

ACT 

1.22 ±0.31 

3.9 

3.2 (9) 

0.96 

ACTPol 

0.85 ±0.36 

2.4 

7.2 (9) 

0.62 

Comb. 

1.06 ±0.24 

4.5 

11.0 (19) 

0.92 


Table 2. Results showing the bias amplitude A relative to the 
fiducial model of Fig. 3. We also quote the signal-to-noise ratio 
S/N, Chi-squared values at the best-fit x^, the number of degrees 
of freedom u and the probability to exceed this x^ (PTE) under 
the assumption of the best-fitting model. 


reported in Table 1. We find A act = 1.22 ± 0.31 and 
AactpoI = 0.85 ± 0.36, with combined constraint A — 
1.06 ± 0.24. The goodness-of-fit statistics for these best-fit 
models are reported in Table 2. This amplitude is consis¬ 
tent with the expected bias from the radio simulations; we 
interpret the result further in Section 3.1. 

The parameter A only scales the bias-dependent part 
of the theoretical model. To assess the overall detection sig¬ 
nificance we rescale the amplitude of the total theoretical 
cross-spectrum by a free parameter a: Cf® —>■ aCi^. This 
is equivalent to equally rescaling both terms in Eq. 4, in¬ 
cluding the magnification bias term. The combined data re¬ 
quire a = 1.06 ± 0.24, and the cross-spectrum is detected 
y /Xnuii “ Xbf = 4.4(7 statistical significance. Here Xnuii — 
31.4 is the chi-squared value of the fit under the null hypoth¬ 
esis (no cross-correlation) and Xbf = H-O is the chi-squared 
value for the best-fit model (number of degrees of freedom 
V = 19). 

The mean cross-spectrum of ACT and ACTPol Monte- 
Carlo simulations with the FIRST dataset is shown in Fig. 2 
and is consistent with null (Section 2.3). These simulations 
reproduce the amplitude and statistics of the leasing field 
but not the true mass distribution on the sky. 

We further test our pipeline, checking for spurious cor¬ 
relations present only in the leasing and galaxy data, by 
performing two additional null tests. 

First, we randomly permute the six FIRST patches 


Figure 6. Cross-correlation between shuffled FIRST maps with 
ACT lensing convergence. Fitting the normalization of the fiducial 
bias model A to these data, we obtain Aact, shuffle = —0.18±0.31, 
consistent with null (Section 3). Grey solid curve: Cross-spectrum 
for the fiducial bias model which best fits the data of Fig. 4. 

within the equatorial strip, such that all patches are moved 
from their true position, with respect to the fixed ACT 
patches. We recompute the cross-spectrum, shown in Fig. 6. 
Fitting the normalization of the fiducial bias model A to 
these data, we obtain Aact, shuffle = —0.18 ± 0.31, consis¬ 
tent with null. The chi-squared value of the null hypothesis 
is 18.2 for v — 10 degrees of freedom, or a probability-to- 
exceed the observed chi-squared of 5%. 

Second, we make reconstructions of the lensing field 
where the deflection field has been redefined - as the curl 
of the lensing potential - and hence the expected ‘conver¬ 
gence’ Q is zero, following Sherwin et al. (2012) and van En- 
gelen et al. (2014). These maps contain reconstruction noise 
but should contain no common signal with the overlapping 
galaxy field. We recompute the cross-spectrum of the lens¬ 
ing curl maps Q with the FIRST maps, shown in Fig. 7. Fit¬ 
ting the normalization of the fiducial bias model A to these 
data, we obtain An = 0.19 ± 0.17. Error bars are calculated 
from the data auto-spectra using the Knox formula (Knox 
1995). The chi-squared value of the null hypothesis is 21.8 
for u = 19 degrees of freedom, or a probability-to-exceed of 
0.29, confirming a null result. 

Removal of the known z < 0.2 sources, which constitute 
~ 5% of the FIRST sample (Section 2.2), has only a small 
effect on the inferred bias amplitude: without removal we 
find a combined constraint AnoZcut = 1.08 ±0.24, consistent 
with expectations given the shape of the lensing kernel and 
low bias of SFGs at low redshift. 

3.1 AGN Bias 

The fiducial bias model scaled by our best-fit amplitude is 
shown in Figure 8. We determine the redshifts to which 
our measurement is most sensitive by considering the kernel 
C^^(z) = WK(z)Wg(z)P(l/x(z); z) of the theoretical cross¬ 
spectrum (Eq. 3; shown in Fig. 9). At I = 200, where the 
signal-to-noise of the cross-spectrum peaks, the mean red- 
shift of the kernel is f zC 2 oo{z)dz/C 2 oo = 1-5. We adopt 
Zeff = 1.5 as the effective redshift of the measurement, es¬ 
timating 6(zefr) = 3.5 ± 0.8. We note however that we are 
sensitive to a range of redshifts: at I = 200 the kernel is 
non-negligible (> 10% of its peak value) out to redshifts 
z > 3, and the kernel shifts to higher redshifts at smaller 
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Figure 7. Cross-spectrum between ACT-l-ACTPol lensing 
curl maps H and FIRST (Section 3). We restrict to Z < 2000 for 
comparison with Fig. 4. Fitting the normalization of the fiducial 
bias model A to these data, we obtain = 0.19 it 0.17 and 
a chi-squared value of the null hypothesis of 21.8 for u = 19 
degrees of freedom (a probability-to-exceed of 0.29). As expected 
this cross-correlation is consistent with null. Grey solid curve: 
Cross-spectrum for the fiducial bias model which best fits the 
data of Fig. 4. 



Figure 9. Cross-power spectrum kernel = 

WK(z)Wg(z)P{l/x(z)-, z), demonstrating the scale-dependent 
sensitivity of the cross-spectrum to source redshift. At I = 200, 
where the signal-to-noise peaks, the mean redshift of the kernel 
is f zC2QQ(z)dz/C20Q = 1.5, which we adopt as the effective 
redshift z^g of the radio source bias measurement. The spread 
in the kernel reflects sensitivity to a wide range of redshift. See 
Section 2.3 for details. 



Figure 8. Solid lines: Fiducial bias model scaled by the best- 
fit bias amplitude for ACT X FIRST (blue), ACTPol x FIRST 
{green) and for the combined result (black). The grey shaded area 
shows the corresponding Icr credible region for the combined re¬ 
sult (Table 3). Note we constrain only the normalization of the 
fiducial bias model; the grey band is indicative of this uncertainty 
only. Using the combined measurements the bias at the effective 
redshift of the analysis is b(zeff = 1-5) = 3.5it0.8 (Section 3). Data 
points: Bias inferred from sources with known individual redshifts 
(left) and with unknown individual redshifts (right), plotted at the 
effective redshift of each measurement (Section 3.1). 

scales. High-redshift sources make an important contribu¬ 
tion to the small-scale cross-spectrum. 

We consider a set of variations to the bias model. First 
we fit the data with a redshift-independent bias model, vary¬ 
ing the amplitude b. We find b = 2.8 ± 0.6, with Xhf = 11-2. 
Our data alone cannot distinguish a redshift-independent 
bias model from an evolving bias model, although our 
redshift-dependent model is more physically motivated by 
theoretical and empirical observations (e.g., Wilman et al. 
2008; Lindsay, Jarvis & McAlpine 2014). 

We probe the flux dependence of the AGN bias by 


splitting the FIRST sources into two roughly equal-sized 
subsamples, with a partition at 2.5 mjy. We create new 
maps of these FIRST sources, as described in Section 2.2. 
From SKADS we estimate the normalised redshift distribu¬ 
tion for each subsample, finding they are equal to within 
~ 15% across 0.3 < z < 4, calculate Zeg = 1.5 for both 
subsamples, and estimate the bias amplitude is ~ 20% 
higher in the bright sample. We repeat the analysis of Sec¬ 
tion 2.3, hnding b(zeg-F > 2.5 mJy) = 4.0 ± 1.1 and 
b{zeH',F < 2.5 mJy) = 3.0 ± 1.1. This is consistent with 
the expectation that the high-flux sample selects preferen¬ 
tially for the most-luminous sources, and these sources lie 
in the most highly-biased environments (e.g.. Mo & White 
1996; Raccanelli et al. 2012). 

We investigate whether the data provide information 
about the (largely unknown) high-redshift bias of radio-loud 
AGN. Here we fix the bias to the fiducial model at redshifts 
z < 1.5 and to a redshift-independent value above this. We 
constrain the high redshift bias to be b{z > 1.5) = 4.1 ± 1.7. 
The detection significance is reduced relative to the full sam¬ 
ple because only high-redshift sources (« 1/3 of the total) 
constrain this parameter. The increase in bias between low 
and high redshift samples is consistent with the result of 
Lindsay, Jarvis & McAlpine (2014), who show that the bias 
b{z) continues to increase above redshift z = 2, although we 
note the significance is low. 

We divide the source sample into those that have red¬ 
shift estimates or not. A fraction of 0.27 of the FIRST 
sources have a reliable optical match as described in Sec¬ 
tion 2.2. The redshift distribution of these sources is strongly 
weighted to low redshifts, peaking around z = 0.5. Follow¬ 
ing Lindsay, Jarvis & McAlpine (2014) we can estimate the 
redshift distribution of the remaining sources by compari¬ 
son with the SKADS simulated radio catalogue used to con¬ 
struct the model redshift distribution for the full sample. We 
construct independent overdensity maps for these two ra¬ 
dio populations (with/without redshift) and recompute the 
data cross-spectra. We also recalculate the theoretical cross¬ 
correlation curves, accounting for the different source distri¬ 
butions, as a function of a redshift-independent bias term b. 
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We find b — 2.1 ± 1.1 at an effective redshift z^b = 0.5 for 
the sample with redshifts, and fe = 3.1 ± 0.8 at an effective 
redshift Zes = 1.6 for the sample without redshifts, shown in 
Fig. 8 . Although not formally significantly different, this is 
consistent with an increasing bias as a function of redshift. 

3.2 Comparison to previous bias measurements 

Geach et al. (2013) find a constant linear bias b = 1.61±0.22 
at an effective redshift z ~ 1 for IR-selected quasars from 
WISE in cross-correlation with the SPT convergence map. 
At the same effective redshift our bias amplitude determina¬ 
tion corresponds to b{z = 1) = 2.6 ± 0.6. Their quasar sam¬ 
ple is shallower (42 sources deg“^) than in the FIRST maps 
presented here (71 sources deg“^), and the predominant sig¬ 
nal comes from z < 2 sources (there are expected to be no 
z > 3 sources). The higher bias determination presented 
here is consistent with a more highly biased population of 
sources being sampled. 

Sherwin et al. (2012) constrain a constant linear bias 
b = 2.5 ± 0.6 for optically-selected quasars from SDSS to¬ 
talling 75 sources deg“^. The redshift distribution of these 
sources peaks at z = 1.4. White et al. (2012) determine b — 
3.8 ± 0.3 from the two-point correlation function of quasars 
in the Baryon Oscillation Spectroscopic Survey across the 
redshift range 2.2 < a < 2.8. Comparing with Fig. 8 , this is 
in good agreement with our result and assumed bias model. 

Lindsay, Jarvis & McAlpine (2014) measure the bias as 
a function of redshift by auto-correlation of radio sources 
from the GAMA survey to the same depth (1 mjy) as this 
analysis. Assuming comoving clustering, their low redshift 
measurement, b{z ~ 0.5) = 2.13to'7g, is consistent with 
the results presented here, while at high redshift they find 
b{z ~ 1.5) = 9.45lg g 7 , significantly higher than seen in this 
analysis. 

Our result probes the multipole range 100 < I < 3000, 
corresponding to physical scales « 2-60 Mpc at the effective 
redshift Zes = 1.5. As seen in Fig. 4, low- and high-redshift 
sources contribute to the cross-spectrum differently as a 
function of scale. At the scales probed by the Planck Collab¬ 
oration (2014b) lensing cross-correlation analysis with NVSS 
radio sources {I < 400), a > 1.5 sources contribute ~ 1/3 
of the total cross-spectrum, whereas at smaller scales these 
sources contribute equally alongside the 2 < 1.5 sources. By 
measuring the cross-spectrum across a wide-range of scales 
one might distinguish between low and high redshift sources. 
Future high-precision determinations of this cross-spectrum 
will go further in breaking the degeneracy between source 
populations and constraining the bias as a function of red¬ 
shift. 

We can translate the constraint on the AGN bias at 
redshift 2 efi = 1.5 into an inference on the mass of the halo 
in which the typical AGN source resides. Using the fitting 
function of Tinker et al. (2010), we find log(Mhaio/Af 0 ) = 
13.6/Io' 4, assuming that haloes virialize at a density ratio 
A = 200 times that of the universe at the epoch of forma¬ 
tion. This observed mass is higher than seen in e.g., Sherwin 
et al. (2012): log(Mhaio/Af 0 ) = 12 . 9 lg; 5 ; Geach et al. 
(2013): log(Mhaio/(/i~^M 0 )) = 12 . 3 lg; 2 . This is consistent 
with the observations of e.g., Shen et al. (2009) andHatch 
et al. (2014) that the environments of radio-loud AGN are 
signihcantly denser than for radio-quiet AGN. 


We hnd a high bias for these sources compared to 
optically- and IR-selected AGN. This analysis provides com¬ 
plementary information by probing the bias of radio-selected 
AGN which, in the context of previous work, is indicative of 
bias evolution and a very large halo mass for these sources. 
The broad picture is that of an increasing bias as a function 
of redshift, and of radio-loud AGN occupying more massive 
haloes than radio-quiet AGN across a similar redshift range. 

Our Endings are in line with studies of the stellar masses 
(e.g., Jarvis et al. 2001b; Seymour et al. 2007) and environ¬ 
ments (e.g., Wylezalek et al. 2013; Hatch et al. 2014) of 
powerful radio sources to high redshift. Specifically, we find 
strong evidence that powerful radio sources are more highly 
biased tracers of the dark matter density field than other 
AGN that are detectable to high redshift (e.g., quasars; 
Sherwin et al. 2012; Geach et al. 2013). As well as being 
important for tracing the underlying dark matter distribu¬ 
tion with techniques such as described in Ferramacho et al. 
(2014), this also suggests that mechanical feedback from the 
jets of powerful radio AGN, should only have a significant 
effect on the level of star formation within the most massive 
dark matter haloes at all epochs. However, we note that such 
an effect can not only have an impact on both the AGN host 
galaxy (e.g., Croton et al. 2006; Bower et al. 2006; Hopkins 
et al. 2006; Dubois et al. 2013; Mocz, Fabian & Blundell 
2013), but also the wider cluster environment (e.g., Rawl¬ 
ings & Jarvis 2004). 

3.3 Modelling limitations and astrophysical 
systematics 

The SKADS simulation is populated using empirical radio 
luminosity functions as described in Wilman et al. (2008). 
Extrapolation of the empirical luminosity functions into un¬ 
observed regimes will lead to inaccuracies in the inferred red¬ 
shift distribution and bias model. To investigate the sensi¬ 
tivity of our measurement to uncertainties about the source 
redshift distribution, we recompute the theoretical spectra, 
unrealistically removing all sources above redshift 2 > 3 
when calculating dn/dz\ at high redshift the underlying 
dn/dz is most uncertain and likely depends on radio lumi¬ 
nosity (e.g., Jarvis & Rawlings 2000; Jarvis et al. 2001b; Wall 
et al. 2005; Rigby et al. 2011). Fitting the theoretical cross¬ 
spectrum as in Section 3, we find b{zeB = 1.2) = 3.2 ± 0.8, 
representing a ~ 0.25a shift from the primary result un¬ 
der this significant perturbation of the theoretical redshift 
distribution. We thus do not expect that the source distribu¬ 
tion uncertainty strongly biases our result, although future 
analyses with higher statistical power will require careful 
consideration of this systematic uncertainty. 

We fix the cosmology to the Planck best-fit values 
throughout this analysis, which could affect the inference 
of the AGN bias. However, the significant (40 ct) detection 
of the Planck lensing auto-spectrum means that model un¬ 
certainty from the cosmology is sub-dominant with respect 
to astrophysical uncertainties (Planck Collaboration 2015). 
Perturbing the Planck best-fit cosmological parameters by 
-l-lcr and recomputing the theoretical cross-spectrum, G/®, 
the amplitude is shifted by < 6% across all relevant scales; 
we thus neglect this source of systematic uncertainty. 

Potential astrophysical systematic contaminants in¬ 
clude infrared sources, Sunyaev-Zeldovich (SZ) clusters and 
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Galactic cirrus. Sherwin et al. (2012) show that these con¬ 
stitute small effects on the measured cross-spectrum be¬ 
tween quasars and lensing (< 10% in total), negligible at 
the level of statistical uncertainty in this analysis. Although 
the sources studied in Sherwin et al. (2012) are optically- 
selected AGN, we expect the result to hold for the radio-loud 
AGN of this analysis given the roughly similar redshift dis¬ 
tributions. Furthermore, bright radio sources (> 5 mjy) in 
the CMB temperature and polarization maps are subtracted 
prior to lensing reconstruction, using a match-filtered source 
template map, thus mitigating radio-source contamination 
in the CMB convergence map (Das et al. 2011; van Engelen 
et al. 2014). 


4 CONCLUSIONS 

We present a measurement of the angular cross-power spec¬ 
trum between lensing convergence from ACT and the over¬ 
density of radio sources identified in the FIRST survey, re¬ 
jecting the null-hypothesis of no correlation at 4.4 ct signif¬ 
icance. The data are well fit by the Planck best-fit ACDM 
cosmological model where we model the source population 
with a redshift-dependent bias. We interpret the result in 
terms of a constraint on the bias of AGN, which dominate 
the FIRST sample, considering various bias models and data 
splits to probe different redshift regimes and AGN popula¬ 
tions, and put these in the context of previous measurements 
of AGN bias. We translate the bias determination into a con¬ 
straint on the mass of the host haloes, corroborating previ¬ 
ous work showing that the environments of radio-loud AGN 
are more dense than those of optically-selected AGN. 

We consider various sources of systematic uncertainty, 
both astrophysical contaminants and modelling limitations. 
We conclude that our results are robust to these effects. 
As deeper and wider radio surveys and improved lensing 
maps become available, these systematic effects will become 
increasingly important to measure and model accurately. 
The auto- and cross-spectra {Cf^,CCj provide com¬ 

plementary information about the large-scale structure they 
probe, with the cross-spectrum in particular being robust 
to systematic biases particular to each dataset. A full anal¬ 
ysis will simultaneously estimate the three power spectra, 
marginalizing over uncertainty in the redshift distribution 
and cosmology (Pearson & Zahn 2014). With current data 
there are strong degeneracies in the cross-spectrum ampli¬ 
tude between sources from different redshifts. The shape of 
the power-spectra contain information about the bias evo¬ 
lution, and larger, more sensitive surveys will allow us to 
break these degeneracies. 

The measurement of the high bias (and correspond¬ 
ingly large halo mass) of this radio population, relative 
to other dark matter tracers, indicates that these sources 
would be useful in the multi-tracer technique of Ferrama- 
cho et al. (2014). Using all the information in auto- and 
cross-correlations between multiple tracers, which differen¬ 
tially trace the dark matter, will provide tight constraints 
on primordial non-Gaussianity by reducing the impact of 
cosmic variance at large scales. 

The SKA will serve as a deep probe of large-scale struc¬ 
ture in the universe, it will be limited by different systemat- 
ics than optical surveys, and the observed source distribu¬ 


tion will be skewed to higher redshifts than either LSST or 
Euclid (Jarvis et al. 2015). Kirk et al. (2015) show that next- 
generation GMB lensing experiments, in combination with 
the SKA, will constrain the amplitude of the lensing-radio 
density cross-spectrum to the sub-percent level. With tight 
constraints on cosmology, this translates into < 1% uncer¬ 
tainty on the bias amplitude, offering broad scope for prob¬ 
ing the history and evolution of AGN. Future high-precision 
measurements of Gf® will use information about the shape 
of the cross-spectrum, and source tomography, to constrain 
the bias as a function of redshift, calibrating galaxy redshift 
surveys and constraining extensions to AGDM. 

The cross-correlation of CMB lensing with tracers of 
large-scale structure will become an increasingly important 
calibrator for future high-precision galaxy and weak-lensing 
surveys. 
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